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ABSTRACT 

The aim of this study is to investigate systematic chemical differentiation of molecules 
in regions of high mass star formation. We observed five prominent sites of high mass 
star formation in HCN, HNC, HCO+, their isotopes, C^*0, C^'^S and some other 
molecular lines, for some sources both at 3 and 1.3 mm and in continuum at 1.3 mm. 
Taking into account earlier obtained data for N2H+ we derive molecular abundances 
and physical parameters of the sources (mass, density, ionization fraction, etc.). The 
kinetic temperature is estimated from CH3C2H observations. Then we analyze corre- 
lations between molecular abundances and physical parameters and discuss chemical 
models applicable to these species. The typical physical parameters for the sources in 
our sample are the following: kinetic temperature in the range ~ 30 — 50 K (it is 
systematically higher than that obtained from ammonia observations and is rather 
close to dust temperature), masses from tens to hundreds solar masses, gas densities 
~ 10^ cm~'^, ionization fraction ^ 10^^. In most cases the ionization fraction slightly 
(a few times) increases towards the embedded YSOs. The observed clumps are close to 
gravitational equilibrium. There are systematic differences in distributions of various 
molecules. The abundances of CO, CS and HCN are more or less constant. There is no 
sign of CO and/or CS depletion as in cold cores. At the same time the abundances of 
HCO"*", HNC and especially N2H"'" strongly vary in these objects. They anti-correlate 
with the ionization fraction and as a result decrease towards the embedded YSOs. For 
N2H''" this can be explained by dissociative recombination to be the dominant destroy- 
ing process. N2H+, HCO"*^, and HNC are valuable indicators of massive protostars. 
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1 INTRODUCTION 

It is now well established that the central parts of dense 
low mass cloud cores suffer strong depletion of molecules 
onto dust grains. Best studied is CO, which has b een shown 
to be depleted in, e. g., L1544 ("Casel li et aLlll999l), IC 5146 
llKramer et all Il999l ). L 1498 (Willa cy. Langer &: Velusanivl 
Il998l ). and L1689B (|jessop fc Ward-ThompsonI I2OOII ). 
Species related to CO, such as HCO^, are als o expected to 



disap pear at gas d ensities above ~ 10 cm ( Caselli et al 



Bergin et al 



Moreover, iTafalla et all (|2002l ') and 
have shown that CS also depletes out in the central 
parts of dense cores, suggesting that CS (so far considered a 
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high density tracer) does not actually probe the central core 
regions. On the other ha nd, N2H^ is an exc ellent tracer of 
dust continuum emission l|Caselli et al.|[200^ ). implying that 
this species does not deplete out (due to the volatility of the 
parent species N2). 

In more massive cores, depletion is probably active in 
dense regions away from star forming sites where dust tem- 
peratures may be low enough (T < 20 K) for CO and 
CS abundances to d rop and cause chemical differentiation 
(|Fontani et al.ll2006l ). 

Several years ago we mapped several tens 
of dense cores towards water masers in CS(2- 
1) with the SEST-15m and Onsala 2 0m radio 

telescopes (IZinchenko. Mattila fc Toriseval Il995l : 

IZinchenko. Pirogov fc Toriseval 119981 ). In 2000 many of 
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them were mapped in N2H+ (|Pirogov et all |2003| ). The 
goal was to identify dense clumps as local maxima in N2H"'" 
maps and to further investigate their properties. However, 
large differences between the N2H"'" and CS distributions 
have been found. 

In this situation it is important to understand which 
species better trace the total gas distribution: is it N2H''' 
as in low mass cores? What is the reason for this differ- 
entiation in warm clouds, where freeze-out is hardly effec- 
tive? To answer these questions we observed dust continuum 
emission and several additional molecular lines towards se- 
lected sources which show significant differences between the 
CS and N2H^ maps. The results for the southern sources 
(where we observed only N2H"'" J — 1 ~0, CS J — 2—1 and 
J = 5 — 4 and dust contin uum) have been published sepa- 
rately (jPirogov et al.l [20071 . hereafter Paper I). In that pa- 
per we have shown that the differences in the CS and N2H'*' 
maps cannot be explained by molecular excitation and/or 
line opacity effects but are caused by chemical differentia- 
tion of these species. We found that N2ll^ abundance in 
many cases drops significantly towards embedded luminous 
YSOs. However, the reasons for this behavior were not clear. 
Two possible explanations were mentioned: an accelerated 
collapse model suggested bv lLintott et al.l (|2005t ) and disso- 
ciative recombination of N2H"'". 

Here we present and discuss the results for the northern 
sample where we observed also HCN, HNC, HCO^, their 
isotopes, C^^O and some other molecular lines, for some 
sources both at 3 and 1.3 mm. These data help to under- 
stand better the chemical differentiation in these objects. In 
addition they give important information on their physical 
properties. 



2 OBSERVATIONS 
2.1 Sources 

The sources for this investigation were selected from 
the sample of m assive cores studied by us earlier 
in various lines ([Zi nchenko. Pirogov fc Toriseva| 1 19981 : 
[Zinchenko. Henkel fc M ao 2000; Pirogov ct al. 20o5). The 
main criterium for this selection was a presence of significant 
differences in the CS and N2H^ maps. The list of the sources 
is given in Table[T] For S187, W3 and S140, the source coor- 
dinates correspond to water masers. For S255, the IRAS po- 
sition is used. In the DR-21 area we used the coordinates o f 
an ammonia core as given bv [jiiina. Mvers fc Adamj ([l999l ') 
which practically coincide (within a few arcsec) with the po- 
sition of DR 21 (OH). In the last column the distances to 
the sources are indicated. The data sets obtained for these 
sources are different. 



2.2 Instruments and frequencies 

The sources were observed with the 20-m Onsala, 12-m 
NRAO (which belongs now to the Arizona Radio Observa- 
tory) and 30-m IRAM radio telescopes in the 3 mm and 1.3 
mm wavebands. Several molecular transitions were observed 
in each waveband. At IRAM 30m also the dust continuum 
emission was mapped at 1.2 mm. The details of observa- 
tions at each instrument are given below. A part of these 



Table 1. Source list. 



Source 


a(2000) 


(5(2000) 


D 




m s J 


(° ' ") 


(kpc) 


S 187 


01 23 15.0 


61 48 47 


1.0 " 


W3 


02 25 28.2 


62 06 58 


2.1 * 


S 255 


06 12 53.3 


17 59 22 


2.5 


DR-21 NH3 


20 39 00.4 


42 22 53 


3.0 


S 140 


22 19 18.2 


63 18 49 


0.9 



Fich fc Blit j lll984h. llBritz Fich fc StarkI l|l982h . 
Harvey. Campbell fc HofFmanrJ lll977t) 

observations has been alre ady published (s o me o f the CS 
and C^'^S J = 2 - 1 data bvlZinchenko et al.J ([l998l ) and the 
N2H+ J = 1 - data by Pirogov ot al. (200J)). We express 
the results of the line observations in units of main beam 
brightness temperature (Tmb) assuming the main beam effi- 
ciencies (77mb) as provided by the telescope documentation. 

2.2.1 Onsala observations 

The observations were performed with SIS receiver in a 
single-sideband (SSB) mode using either dual beam switch- 
ing with a beam throw of 11'5 or frequency switching. As a 
backend until the year 2000 we used 2 filter spectrometers 
(usually in parallel): a 256 channel filterbank with 250 kHz 
resolution and a 512 channel filterbank with 1 MHz reso- 
lution. Since 2000 we used mainly the autocorrelator spec- 
trometer tuned to 50 kHz resolution. Pointing was checked 
periodically by observations of nearby SiO masers; the point- 
ing accuracy was typically < 5". The half-power beam width 
(HPBW) was from about 35" at the highest frequencies to 
about 40" at the lowest frequencies. 

The standard chopper-wheel technique was used for the 
calibration. The system temperature varied in a wide range 
depending on the weather and observing frequency, from 
~ 200 K at lower frequencies in good weather to ~ 1000 K 
and more at higher frequencies and under cloudy conditions. 

2.2.2 NRAO 12m observations 

At the NRAO 12-m telescope only two sources from those 
listed in Table[T]were observed: S187 and S255 (in the C^**0, 
CS, C^'^S, SiO and methanol lines at 1.3 mm). As a result, 
the data sets for these sources are the most complete ones. 

The observations were perfomed in 2000 with the SIS 
receiver and two backends in parallel: the filter bank with 
the 0.5 MHz resolution and MAC autocorrelator with the 
100 kHz resolution. We used frequency switching and po- 
sition switching observing modes. The pointing and focus 
were checked periodically by observations of planets. The 
HPBW was 26-27". The system temperature varied from 
~ 300 K in good weather conditions to ~ 1500 K with in- 
creasing humidity. 

2.2.3 IRAM 30m observations 

At the IRAM 30m telescope we obtained maps of the con- 
tinuum emission from our sources at 1.2 mm and their maps 
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Table 2. Results of the continuum observations: brightness (B), 
fluxes (F) and angular sizes (9) of the emission clumps towards 
selected positions. 



Source 


Aa 


AS 


B 


F 


e 


Peak 




{") 


{") 


(Jy/beam) 


(Jy) 


(") 




S 187 


+160 





0.15 


2.1 


45 


CS 







+80 


0.09 


0.46 


28 


N2H+ 


W3 


20 


-40 


2.03 


19 


38 


CS 




+160 


-160 


0.70 


2.23 


22 


N2H+ 


S 255 








1.35 


6.4 


27 


CS 







+60 


1.24 


6.4 


28 


N2H+ 


S 140 








1.58 


15.9 


39 


CS 




+40 


+20 


0.34 


5.3 


49 


N2H+ 



in several components of the CH3C2H J = 13 — 12 transi- 
tion. The continuum observations were performed with the 
MAMBO bolometer array (the details of this instrument are 
available at the IRAM Web site) and reduced with the MOP- 
SIC package. The spectral line observations were done with 
the HERA (Heterodyne Receiver Array) and VESPA auto- 
correlator backend. The pointing and focus were checked pe- 
riodically on nearby strong continuum sources. The antenna 
HPBW is about 12" at these frequencies. The typical sys- 
tem temperatures for HERA observations were in the range 
~ 200 - 400 K. 



3 OBSERVATIONAL RESULTS 

We present the observational results in the form of maps 
and tables where the line parameters at selected positions 
are given. These positions correspond to different emission 
peaks in the sources which are identified in the maps of 
molecular emission and in most cases can be seen in the 
continuum maps presented in Fig. [T] These maps are plotted 
using logarithmic scale for intensity in order to emphasize 
weak features. The continuum brightness, fluxes and angu- 
lar sizes of the emission clumps towards selected positions 
are summarized in Table [2] The sources are rather extended 
in continuum. In order to provide a better comparison with 
the molecular data, most of which were obtained with ~ 30" 
(HPBW) beams, we give the fluxes integrated over 1' circles 
centered at the selected positions. These positions are la- 
belled as "CS" and "N2H^" emissi on peaks according to 
our previous CS an d N2H+ surveys (jZinchenko et al.lll998l : 
IPiroeov et al]l2003h and other available data. The sizes are 
derived from these fluxes and brightness as 6 = ^ AF/ixB. 
For Gaussian brightness distribution this corresponds to size 
at the 1/e level which is about 20% larger than the size at 
the half maximum level (^0.5 )• In S140 it is hard to see a 
distinct continuum clump at the N2H+ emission peak, one 
can see rather an extended filament here. Nevertheless, we 
provide brightness and flux for this position, too. The spe- 
cific features of every source are briefly described below. The 
spectral data reduction was performed with the GILDAS 
and XSpec (Onsala data) software packages. 



3.1 Molecular maps and line parameters 

3.1.1 S187 

Maps of S187 in various lines overlaid on the greyscale map 
of the 1.2 mm continuum emission are presented in Fig. [J] 
The CS(2-1) and N2H"'"(l-0) maps have been pu blished ear- 
lier ()Zinchenko et al.lll998l : IPiroeov et al.ll2003l ). There is a 
striking difference between the various maps. Two separated 
clumps are clearly seen in N2H"''(l-0), which peaks North- 
West of the IRAS 01202+6133 source. The two peaks are 
still visible in the HNC(l-O) map, although they are not 
as well separated as in the N2H^(l-0) map. All the other 
maps have a cometary or more irregular morphology, with 
peaks all shifted from the continuum peak. The gaussian 
line parameters at the CS and N2H+ emission peaks are 
summarized in Table [3] In C^*0 J=l-0andJ = 2- l 
emission towards the CS peak there is an additional narrow 
(~ 0.7 km/s) component at about —15.5 km/s. However, it 
is not pronounced in the lines of other high density tracers 
and we do not include it in Table |3] 

The molecular data indicate the presence of at least 3 
clumps in the area. There are several IRAS point sources and 
molecular masers here. The strongest IRAS point source, 
IRAS 01202+6133, is located at about 2' to the east from 
our central position and coincides with the main 1.2 mm 
continuum p eak. Here an OH maser and UC H 11 region 
are present (| Argon. Reid fc MentenI I2OO0I ). The secondary 
N2H"'" peak coincides with this IRAS position. Also a weaker 
CS clump is located here as is clearly seen from the CS J — 
5 — 4 data. The main CS, HCO^ and HCN emission peaks 
are shifted by about 1' further to the east. No IR sources or 
masers are known in this area. It is worth noting that the 
methanol emission peak is shifted still further to the east. At 
the same time C^**0 emission peaks near IRAS 01202+6133. 

The strongest N2H^ peak coincides with a rela- 
tively weak 1.2 mm continuum clump. It is shifted by 
about 0!5 from the strong near IR source NIRS 60 
(|Salas. Cruz-Gonzalez fc Porraslll998t ). 



3.1.2 W3 

The molecular line maps of W3 are shown in Fig. [3] 
overlaid with the 1.2 mm continuum emission map. The 
CS(2-1) and N2H"'"(l-0) maps have been p ublished earlier 
jZinchenko et al.lll998l : IPirogov et al.ll2003r ). The brightest 
continuum peak is attrib uted to free-free emis sion from a 
compact H 11 region (e.g. iTieftrunk et"aLlll997h . There are 
two main molecular emission peaks in this area associated 
with two dust clumps. As in the S187 case, the N2H^ map 
looks very different from most other maps. The N2H"'' peak 
coincides with a relatively weak south-east (SE) clump at 
the approximately (160",— 160") position. The gaussian line 
parameters at the CS and N2H"'" emission peaks are summa- 
rized in Table|31 Some spectra clearly show non-gaussian fea- 
tures: broad wings at the CS peak position and red-shifted 
self-absorption in HCO"^ at the N2H"'" peak. 

3.1.3 S255 

For S255 we have the most complete data set. In addi- 
tion to the common set of molecular transitions it includes 
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Aa (") Aa (") 

Figure 1. Maps of the sample sources in continuum at 1.2 mm obtained at the IRAM 30-m telescope. The sources are (from left to 
right and from top to bottom): S187, W3, S255 and S140. The intensity scale is logarithmic. The triangles mark water masers, the stars 
indicate point IR sources (IRAS point sources for S18 7, S255 and W3; IRSl, IRS2 and IR S3 for S140). The circle corresponds to the 
submillimetre source SMAl in the S255 N clump area jCyganowski. Brogan &: Huntei]|2007l) . 



SO(32 - 2i), SiO (2-1) and (5-4), several methanol 2-1 
and 5-4 series li nes. Earlier we mapped it also in amm omia 
(1,1) and (2,2) (jZinchenko, Henning fc Schreve j[l997l ') . The 
molecular maps overlaid on the greyscale map of the 1.2 mm 
continuum emission are presented in Fig. |4l This source was 
m apped at 1.2 mm in co ntinuum at the IRAM 30m telescope 
bv lMezger et al.l (|l988l ). Our observations with the new ar- 
ray receiver provide a better sensitivity and a wider map 
area, although the basic features of our map are consistent 
with those previous results. 

The maps show two main peaks of molecular and dust 
continuum emission, around the (0,0) and (0,-1-60") posi- 
tions. The commonly accepted names for these clumps are 
S255IR and S255N. There is also a third (southern) peak 



at about (-|-20",— 50") noticeable in the continuum map 
and in several molecular maps, at least in N2H^(l-0) and 
HCO^(l-O). N2H^ and ammonia are significantly stronger 
at the northern peak (S255N) while most other species are 
either stronger at the central position (S255IR) or compara- 
ble at both central and northern clumps. The dust emission 
is almost equal for these clumps. The nature of these two 
components is different. The central one is associated with 
a luminous cluster of IR sources, whereas toward the north- 
ern one an ultracompact H ll region (G192.58-0.04) was de- 
tected. Recently several compact submil limetre continuum 
clum ps were detected there with SMA (jCveanowski et al] 
[2007). T his object is ex t remely red in the m id-IR band 
(jCrowther fc Contil [20o3 ) . iMezger et all (|l988l ) derived al- 
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Figure 2. Maps of S187 in various molecular lines (contours) overlaid on the map of 1.2 mm dust continuum emission (grayscale). The 
contour levels span from 20% to 100% of the peak intensity in steps of 10%. 



Table 3. Molecular line parameters at the CS and N2H+ emission peaks in S187. 





(+160",0) 






(0,+80") 






Line 


Tmb 


Vlsr 


AV 




Vlsr 


Ay 




(K) 


(km/s) 


(km/s) 


(K) 


(km/s) 


(km/s) 


CiSO(l-O) 


5.47(18) 


-13.92(03) 


1.71(09) 


3.24(18) 


-13.78(04) 


1.67(10) 


Ci80(2-1) 


5.85(19) 


-14.05(02) 


1.69(05) 








CS(2-1) 


4.29(10) 


-14.26(03) 


2.30(07) 


1.20(20) 


-13.63(27) 


3.24(70) 


C34S(2-1) 


1.00(04) 


-14.05(04) 


2.25(11) 








CS(5-4) 


3.58(06) 


-14.09(01) 


1.86(03) 








C34s(5-4) 


<0.4 












HCN(l-O) 


5.14(07) 


-14.31(02) 


2.33(04) 


1.92(08) 


-13.81(04) 


1.63(07) 


Hi3cN(l-0) 


0.50(05) 


-14.17(10) 


1.83(21) 








HCO+(1-0) 


5.54(08) 


-14.49(02) 


2.22(04) 


2.41(08) 


-14.03(03) 


1.85(08) 


Hi3CO+(l-0) 


0.58(07) 


-14.15(09) 


1.56(22) 


1.28(10) 


-13.51(03) 


0.83(07) 


HNC(l-O) 


5.44(12) 


-14.29(02) 


2.27(06) 


5.70(15) 


-13.58(02) 


1.39(04) 


HNi3C(l-0) 


0.33(04) 


-13.88(09) 


1.66(20) 


0.63(05) 


-13.35(03) 


0.76(07) 


N2H+(l-0) 


0.61(08) 


-14.09(07) 


1.45(20) 


1.59(12) 


-13.33(03) 


0.90(06) 
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Figure 3. Maps of W3 in various molecular lines (contours) overlaid on the map of 1.2 mm continuum emission (grayscale). The contour 
levels span from 20% to 100% of the peak intensity in steps of 10%. 



Table 4. Molecular line parameters at the CS and N2H+ emission peaks in W3. 





(0,-40") 






(+160",- 


^160") 




Line 




Vlsr 


AV 


Ti-nb 


Vlsr, 


AV 




(K) 


(km/s) 


(km/s) 


(K) 


(km/s) 


(km/s) 


Ci**O(l-0) 


2.05(20) 


-42.70(27) 


5.32(71) 


2.36(37) 


-38.58(11) 


1.49(27) 


CS(2-1) 


8.82(10) 


-43.09(03) 


4.76(06) 


3.19(03) 


-38.66(19) 


4.43(53) 


HCN(l-O) 


7.96(11) 


-43.08(02) 


3.41(05) 


1.65(13) 


-38.55(40) 


8.08(38) 


H"CN(l-0) 


1.51(07) 


-42.75(12) 


4.62(20) 








HCO+(1-0) 


16.26(07) 


-43.40(01) 


3.86(02) 


4.05(06) 


-38.97(04) 


4.73(09) 


H13CO+(1-0) 


0.88(17) 


-43.11(21) 


2.44(51) 


0.95(18) 


-38.33(19) 


2.02(45) 


HNC(l-O) 


5.36(11) 


-42.92(04) 


3.83(09) 


4.41(12) 


-38.74(04) 


3.12(10) 


HNi3c(l-0) 


0.37(03) 


-42.63(09) 


2.80(21) 


0.17(01) 


-38.38(05) 


2.95(23) 


N2H+(l-0) 


0.66(09) 


-42.22(11) 


2.04(30) 


1.37(10) 


-38.76(10) 


1.75(14) 



most exactly the same dust masses and temperatures for 
both components from their 1.2 mm and 350/im observa- 
tions. The gaussian line parameters at the central and north- 
ern emission peaks are summarized in Table [3 



3.14 DR-21 NH3 

This is another well known site of active star formation. 
The CS, HCN and HCO^ spectra here suffer from a very 
strong red-shifted self-absorption (Fig. [5]), which makes an 
analysis of corresponding maps almost senseless. For this 
reason we present here and discuss only rarer isotopologue 
and N2H"'"(l-0) maps (Fig. (6}. We did not map the dust 
continuum emission in this area, given that it ha s been 
already observed by IChandler. Gear fc Chinil (Il993l ). It is 
easy to see that the N2H"'" and HN^'^C peaks are shifted 
by about 0'5 to the south from the emission peaks of 
other species. The latter peaks practically coincide with the 
main dust emission peak DR21(OH)M in the n otation intro- 
duced by iMangum. Wootten fc Mundvl l|l99ll ). The N2H+ 



and HN^^C peaks lie near a weaker dust peak DR21(0H)S. 
The gaussian line parameters are given in Table [6] 



3.1.5 SI40 

S140 is one of the best studied sites of active star forma- 
tion. In particular it was observed at various i nstruments in 
the H CN, HCO+ and CO isotopic lines (e.g. iPark fc Minhl 
Il995h . Nevertheless, we present here our own data on these 
lines too, to better compare with other molecular maps (all 
maps have been obtain ed at Onsala) . T he N2H^ results have 
been published earlier (|Pirogov et al. II2OO3). 

The maps (Fig. clearly show that HCN and HCO+ 
emissions peak near the (0,0) position while the N2H^ peak 
is shif ted to about (-|- 50",-|-20"). A multi-transitional CS 
study (jZhou et al.lll994l ) shows that CS emission peak is near 
the (0,0) position. HNC is an intermediate case and there is a 
secondary HNC peak near the N2H+ peak. These differences 
in emission distributions are not caused by opacity effects 
because the maps in the lines of rarer isotopic modifications 
of these molecules show the same features (the N2H^ opti- 
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Figure 4. Maps of S255 in various molecular lines (contours) overlaid on the map of 1.2 mm dust continuum emission (grayscale). The 
contour levels span from 20% to 100% of the peak intensity in steps of 10%. 



cal d epth is rather small, ~ 0.4, as shown bv lPirogov et al] 
l2003f ). The gaussian line parameters are presented in Ta- 
ble [7] 



4 PHYSICAL PARAMETERS OF THE 
SOURCES AND MOLECULAR 
ABUNDANCES 

4.1 Kinetic temperatures 

The kinetic temperatures of the sources have been esti- 
mated from t he C H3C2H observations. As shown e.g. by 
iBergin et al.l (|l994l ) this symmetric-top molecule is a good 
"thermometer" for gas densities n > 10* cm~^. Our main 
is to compare the temperatures at the peaks of the 



CS and N2H"'" emission. At first, we derived the tempera- 
tures from the CH3C2H J — 6 — 5 data obtained at Onsala. 
These measurements were done only at the peak positions. 
Later the CH3C2H J = 13 — 12 emission was mapped at 
IRAM with the multibeam receiver. The det ails of the tem- 
perat ure estimates are published elsewhere (|Malafeev et al.l 
|2005| . Zinchenko et al., in preparation). Here we present the 
main results. In S187 the CH3C2H emission was too weak. 

The kinetic temperatures at the CS and N2H^ peaks, 
derived from the Onsala and IRAM CH3C2H data, are listed 
in Table [S] It is easy to see that in most cases there is no 
significant temperature difference between the CS and N2H"'" 
emission peaks, although on average the N2H"'" peaks are 
somewhat colder than the CS ones. 

CH3C2H maps of sufficiently high quality were obtained 
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Figure 6. Maps of DR-21 NH3 in various molecular lines. The contour levels span from 20% to 100% of the peak intensity in steps of 
10%. The cross marks the position of DR21 (OH). 
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Table 5. Molecular line parameters at the CS (central) and N2H+ (northern) emission peaks in S255. 





(0,0) 






(0,+dO ) 






Line 


^ mb 




AT/ 


-'mb 




AT/ 






(km/s) 


(km/s) 




(km/s) 


(km/s) 


Ci8O(l-0) 


1.99(12) 


6.84(17) 


4.31(43) 


1.67(13) 


8.74(14) 


3.71(34) 


Ci«0(2-1) 


6.45(08) 


7.18(01) 


2.88(03) 


5.89(08) 


8.82(01) 


2.95(03) 


CS(2-1) 


13.36(09) 


7.46(01) 


2.78(02) 


9.08(08) 


8.40(01) 


3.16(03) 


0343(2-1) 


1.93(11) 


7.46(08) 


2.78(18) 


1.33(10) 


8.56(12) 


3.20(29) 


CS(5-4) 


10.94(12) 


7.49(01) 


2.73(03) 


7.46(12) 


8.55(02) 


3.32(05) 


C34s(5-4) 


1.13(06) 


7.52(05) 


2.33(12) 


0.75(08) 


9.10(13) 


3.57(26) 


HCN(l-O) 


13.79(11) 


7.33(01) 


3.25(03) 


11.52(11) 


8.40(02) 


3.60(03) 


Hi3CN(l-0) 


1.41(07) 


7.57(06) 


2.58(12) 


1.37(07) 


8.91(06) 


2.12(11) 


HCO+(1-0) 


7.34(10) 


7.53(03) 


3.90(06) 


11.27(11) 


8.83(02) 


3.21(04) 


H13CO+(1-0) 


0.74(10) 


7.51(11) 


1.69(26) 


1.39(08) 


8.97(07) 


2.79(18) 


HNC(l-O) 


11.87(18) 


7.30(02) 


2.73(05) 


11.30(17) 


8.55(02) 


3.32(06) 


HN"c(l-0) 


0.18(02) 


7.29(20) 


3.30(47) 


0.28(02) 


8.87(16) 


4.78(37) 


N2H+(l-0) 


1.46(06) 


7.44(04) 


2.25(10) 


2.58(05) 


8.97(02) 


2.61(06) 


SiO(2-l) 


0.18(06) 


7.96(52) 


3.32(122) 


0.20(05) 


8.06(61) 


5.24(144) 



Table 6. Molecular line parameters at the CS and N2H+ emission peaks in DR-21 NH3. 





(0,-20") 






(0,-40") 






Line 




Vlsr 


AV 


2^mb 


Vlsr 


AV 




(K) 


(km/s) 


(km/s) 


(K) 


(km/s) 


(km/s) 


Ci8O(l-0) 


4.99(14) 


-2.99(05) 


3.84(13) 


4.69(15) 


-3.37(05) 


3.28(12) 


C34S(2-1) 


2.62(06) 


-3.04(05) 


4.68(12) 


2.89(06) 


-3.01(04) 


3.83(10) 


Hi3cN(l-0) 


2.34(06) 


-3.24(06) 


4.07(10) 


2.34(06) 


-3.47(05) 


3.58(09) 


Hi3cO+(l-0) 


3.17(07) 


-3.46(04) 


3.67(09) 


2.73(07) 


-3.57(04) 


3.25(10) 


HN"c(l-0) 


1.34(09) 


-3.05(14) 


4.38(32) 


1.70(09) 


-3.25(10) 


3.56(23) 


N2H+(l-0) 


4.99(04) 


-3.11(02) 


4.04(03) 


6.93(04) 


-3.25(01) 


3.51(02) 



Table 8. Kinetic temperatures at the CS and N2H+ peaks de- 
rived from the Onsala and IRAM CH3C2H data (T^j^^ and 
'^khT^^' respectively). 



Source Aa AS T^r^ T^?'^^ Emission 

(") {") (kT (kT peak 



W3 


20 


-40 


52.6 ±3.1 


58.4 ± 1.5 


CS 




160 


-160 


30.7 ±0.8 


37.3 ±2.2 


N2H+ 


S255 





60 


34.9 ± 1.4 


39.6 ±0.4 


N2H+ 










34.5 ± 1.0 


38.8 ±0.9 


CS 


DR21 





-40 


28.8 ±0.1 


34.8 ± 1.1 


N2H+ 










33.4 ±2.7 


44.0 ± 1.4 


CS 


S140 


40 


20 


27.8 ± 1.6 


37.7 ± 2.2 


N2H+ 










30.6 ±0.7 


42.4 ± 1.0 


CS 



It is interesting to compare our estimates of kinetic tem- 
perature with other available data for these sources. Most 
of them have been actively st udied. Such comparison was 
made bv lMalafeev et all (|2005l ) and here we repeat the main 
points. 



4.1.1 W3 

Based on NH3(1,1) and (2,2) lines. iTieftrunk et al.l (|l998l ) 
derived kinetic temperatures of 44 K and 25 K toward our 
CS and N2H+ peaks, respectively. A comparison with Ta- 
ble |5] shows that these temperatures are somewhat lower 
than those obtained from the CH3C2H data, although the 
relation is the same. 



at IRAM for DR21, S140 and the CS peak in W3. The IRAM 
data indicate somewhat higher peak temperatures than ob- 
tained from the Onsala data, as expected due to a higher 
angular resolution and higher excitation requirements for 
the J = 13 — 12 transition if temperature increases towards 
an embedded heating source. Our data do show such tem- 

J ierature grad ients consistent with theoretical expectations 
Zinchenko et a l. 2005) but we do not discuss them here. At 
the same time the IRAM data confirm our main conclusion 
that there is no significant temperature difference between 
the CS and N2H"'" peaks in most cases, although the N2H"'" 
peaks are somewhat colder than the CS ones, on average. 



4.1.2 S255 

From our ammonia observations l|Zinchenko et al.] Il997l ) 
the kinetic temperature at the (0,+80") position (near the 
NH3 and N2H"'" peaks) is 23 ± 1 K. At the central posi- 
tion the uncertainty in the kinetic temperature is too high. 
ISchrever et al.l l|2000 l) obtained the kinetic temperature in 
the centre from their ammonia data of 30.8 K and the dust 
tempe rature from the IRAS data of 31.1 K. iMezger et al.l 
(| 19881 ) derived almost the same dust temperatures (~ 30 K) 
for the two peaks discussed here from their 1.3 mm and 350 
/^m continuum observations. 



© 2009 RAS, MNRAS 000.[TlfT5l 



10 /. Zinchenko, P. Caselli & L. Pirogov 

Table 7. Molecular line parameters at the CS and N2H+ emission peaks in S140. 





(U,UJ 






(+40", +20") 








— — 

Vlsr 






^LSR 


ZA V 




(K) 


(km/s) 


(km/s) 


(K) 


(km/s) 


(km/s) 


Ci8o(l-0) 


3.19(16) 


-6.83(09) 


3.46(22) 


2.53(18) 


-7.51(10) 


2.95(24) 


HCN(l-O) 


18.37(11) 


-6.80(01) 


3.05(02) 


11.43(11) 


-6.92(01) 


2.84(03) 


Hi3cN{l-0) 


1.50(07) 


-7.00(05) 


2.22(10) 


0.79(08) 


-6.96(08) 


1.53(16) 


HCO+(1-0) 


25.63(09) 


-6.91(01) 


3.05(01) 


15.13(09) 


-6.88(01) 


2.86(02) 


Hi3CO+(l-0) 


2.64(10) 


-7.18(05) 


2.58(11) 


1.94(10) 


-7.21(06) 


2.26(14) 


HNC(l-O) 


16.62(17) 


-6.77(01) 


2.83(03) 


14.95(18) 


-7.03(02) 


2.51(04) 


HNi3c{l-0) 


0.68(05) 


-6.42(12) 


3.27(27) 


1.24(07) 


-6.94(05) 


1.75(11) 


N2H+(l-0) 


3.41(07) 


-6.90(02) 


2.20(05) 


6.33(09) 


-7.00(01) 


1.87(03) 




-20 -10 10 

Vlsr (km s""') 



Figure 5. Molecular line spectra towards the N2H+ emission 
peak in DR-21 NH3. 



1.1.3 DR21 NH3 

This source was observed (without mappin g) in the CH3C2H 
lines at the 11-m NRAO telescope (K uiper et al.l Il984l ). 
They observed the same J = 6 — 5 transition and derived 
practically the same kinetic temperature at the centre as in 
our pres ent work, 33.0 K. There are VLA am monia observa- 
tions by[M angum. Wootten fc Mundvl (Il992l ). They derived 
kinetic temperatures of 32 K at the (0,0) position and 25 K 
near the (0,-40") position. This is rather close to our esti- 
mates. 

4.1.4 SI40 

S140 was als o observed in CH3C2H at the 11-m NRAO 
telescope by iKuiper et al.l (|l984l ). They obtained for the 
cloud centre the kinetic temperature of 32.1 ± 6.7 K 
which coincides with our J = 6 — 5 estimate within the 
uncertainties. From ammonia observ a tions at Effelsberg, 
lUngerechts. Walmslev fc Winnewissej (|l986l ) found a tem- 
perature of about 20 K, once again a value significantly lower 
than found w ith CH3C2H. The dust temperature from IRAS 
data is 34 K (|Zhou et al.lll99j ). more similar to our findings. 
Thus, CH3C2H appears to better trace the dust. 

4.2 Masses and densities 

There are several ways to estimate masses and densities of 
interstellar clouds. Here we derive the source masses and 
column densities primarily from the dust continuum obser- 
vations. This is considered to be one of the most reliable 
methods, partly due to the fact that dust/gas mass ratio is 
rather constant. At the same time the uncertainties in the 
dust opacities are still rather high. Additional uncertain- 
ties are related to uncertainties in the dust temperature. 
Here we use the temperature estimates obtained from our 
methyl acetylene observations, found to be similar to the 
dust temperature. It is known that at the typical densities of 
dense cores the dust and gas kinetic temperatures are close 
to each other. The results of these estimates are presented 
in Table |9] It is worth noting that the masses here are not 
really total masses of the sources because they are derived 
from the fiuxes integrated over 1' diameter circle around the 
indicated positions. The dus t opacities were adopted from 
lOssenkopf fc Hennind (|l994l ). The gas/dust mass ratio was 
assumed to be equal to 100. In addition we present also es- 
timates of the average volume density along the line of sight 
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Table 9. Estimates of the gas column densities, average volume 
densities and masses from the continuum observations. In the last 
column estimates of the virial masses are given. 



Source 


Peak 


A^L X 10~2^ 


n X 10"^ 


M 








(cm~^) 


(cm~^) 


(Mq) 


(Mq) 


S 187 


CS 


0.4 


0.7 


30 


50 




N2H+ 


0.5 


1.0 


14 


10 


W3 


CS 


2.0 


1.6 


405 


310 




N2H+ 


1.2 


1.5 


74 


90 


S 255 


CS 


2.0 


1.9 


290 


210 




N2H+ 


1.9 


1.6 


280 


220 


S 140 


CS 


2.8 


3.9 


86 


90 




N2H+ 


0.6 


0.7 


32 


70 



obtained as n = Ni, / L where the size L is derived from the 
angular size (Table [2} and distance to the source (Table [ij. 

In the last column we give estimates of the virial masses 
of the clumps based on these sizes and widths of optically 
thin lines at these positions. These estimates are rather un- 
certain, at least by a factor of 2 due to large uncertainties in 
the parameters. Nevertheless, within this factor they agree 
well with the masses derived from continuum observations. 
This shows that the clumps are close to gravitational equi- 
librium. 

The average gas volume densities are n ~ 10^ cm~^ in 
all cases. However, it is well known that such mean densities 
can be significantly lower tha n densities found from molec- 
ular excitation analysis (e.g. ^Zinc henko et al.lll998l ). This 
discrepancy is most probably explained by small-scale gas 
clur npiness. Such dumpiness is indic ated by many studies 
(e.g. iBergin. SneU fc Goldsmithll 19961 ) . At the same time our 
data do not provide sufficient material for excitation analy- 
sis in most cases. Only for S255 we have 3 mm and 1.3 mm 
data in C^**0, CS, C^^S and methanol transitions towards 
both CS and N2H"'" peaks which can be used for modeling. In 
particular, our LVG estimates of gas density from the C"^'*S 
data give values of about n ~ 5 x 10^ cm~^ for both compo- 
nents. Methanol data modeling gives n ~ 3 x 10^ cm~'^ also 
for both peaks (S. Salii, private communication). Estimates 
based on C^®0 are highly uncertain but are consistent with 
these values and show similar densities for both components 
because the line intensity ratios are similar. These estimates 
are rather close to the average densities presented in Table[51 
although somewhat higher as expected. This shows that the 
volume filling factor in this density range is rather high. 
The continuum data also indicate practically equal mean 
gas densities for the components. 



4.3 Abundances 

We derive here abundances of several species which were 
observed in all our sources and are the most informative 
ones for our purposes: C^*0, H"CN, H"C0+, HN^^'C, 
N2H"'" and C^*S. In most cases we use only rare isotope 
data which are presumably not affected by optical depth 
effects. The optical dept h in the N2H" '" lines is also small 
or moderate as shown bv lPirogov et al.l (|2003D . The column 
densities are estimated from integrated line intensities in 
the LTE approximation assuming the excitation tempera- 



tures equal to the kinetic temperatures as given in Table [8] 
(approximate values close to TJ^j"^ and T^^~^^ were used). 
For S187 we assume Tkin of 20 K for the CS peak and Tkin 
of 10 K for the N2H"*' peak (taking into account narrow 
line widths and very weak CH3C2II emission). Then, we 
derive abundances from comparison of the molecular col- 
umn densities and total gas column densities obtained from 
dust continuum observations (Table [9}. In order to take 
into account different beam sizes we correct the abundance 
estimates assuming Gaussian beams and Gaussian bright- 
ness distributions. For DR21 where we do not have our own 
dust continuum data, the abundances were derived assum- 
ing X(C^ ^O) = 2 X 10~^. Th e results are summarized in 
Table \W\ IPirogov et al] l|2003h derived N2H+ abundances 
in a different way: from the N2H"^ column densities and 
virial masses of the clouds. Nevertheless, their estimates are 
very close to those presented in Table [10] when the positions 
coincide (within a factor of 2). 

Non-LTE modeling for some clumps in our sample using 
the LVG approximation or t he RADEX code (based on the 
escape probability method; Ivan der Tak et ahl |2007^ gives 
systematically lower (by a factor of 1.5-2) column densities 
for all species considered here, except C^*0, assuming tem- 
perature and density as described above. This is a natural 
result because C^^^O is practically thermalized at these con- 
ditions while the other species are still far from equilibrium 
and population of lower levels exceeds that expected in LTE. 
Nevertheless, we believe that non-LTE modeling is not jus- 
tified here due to the absence of the necessary data for most 
of the sources. One can bear in mind that the abundances 
derived here are somewhat overestimated. 



4.4 Ionization fraction 

The ionization fraction can be estimated from the HCO'*' 
abundance using the results of relevant chemical models. It 
is known that in dens e clouds HCO'*' is formed mainly from 
H+ fe.g. lTurneijll995l ): 

-f- CO ^ HCO+ -h H2 (1) 

and is destroyed by dissociative recombination: 

HCO+ + e ^ CO + H (2) 

In addition, one has to take into account HCO^ 
recombination on to negatively charged dust grains (e.g. 
ICaselli et al.ll200^ . The rate of this reaction is given by the 
product fcgXg where the rate coefficient fcg and the fractional 
abundance of dust grains Xg are determin ed by grain size 
distri bution and gas kinetic temperature (IDraine fc SutinI 
[l983). For the typical temperatur es in our objects (~ 
30 K) and MRN l|Mathis et al.lll977l ') grain size distribution 
fcgXg 3 X lO"^'^ s-^ 

In this model for HCO^ abundance in steady state we 
can write: 



X(HCO+) = 



fciX(H+)X(CO) 
a(HCO+)Xe + fcgXg 



(3) 



where k\ is the rate of reaction ([T]) and ci(IICO^) is the 
1100"*" dissociative recombination rate. 

H3" is formed by cosmic-ray ionization of molecular hy- 
drogen a nd is destroy ed in dense clouds primarily by reaction 
g]) (e.g. iBlackllioOol ). For regions where CO is not frozen 
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Table 10. Estimates of molecular and electron abundances towards selected positions in the sample sources. For C^^O, abundances 
derived from the J = 1 — and J = 2 — 1 transitions are shown separately. 



Source Aa A<5 X(Ci**0)io X(Ci80)2i X(Hi3CN) X{Yl^^CO+) X{}m^^C) X{C'-''^S) X{N2H+) Xc 
(") {") 



S187 


160 





4.2E-07 


1.5E- 


-07 


l.OE- 


10 


5.8E- 


11 


6.5E- 


11 


6.0E- 


10 


6.5E-11 


1.7E-07 







80 


2.5E-07 










7.6E- 


11 


6.2E- 


11 






1.2E-10 


1.3E-07 


W3 





-40 


2.9E-07 






4.8E- 


10 


8.6E- 


11 


7.6E- 


11 






6.5E-11 


1.2E-07 




160 


-160 


2.0E-07 










1.9E- 


10 


8.8E- 


11 






2.7E-10 


5.4E-08 


S255 








2.4E-07 


1.4E- 


-07 


2.8E- 


10 


5.7E- 


11 


4.9E- 


11 


9.2E- 


10 


1.7E-10 


1.8E-07 







60 


1.8E-07 


1.3E- 


-07 


2.3E- 


10 


1.9E- 


10 


1.2E- 


10 


7.7E- 


10 


3.7E-10 


5.4E-08 


DR21 





-20 








2.2E- 


10 


1.6E- 


10 


1.5E- 


10 


6.3E- 


10 


3.2E-10 


6.3E-08 







-40 








2.4E- 


10 


1.5E- 


10 


1.9E- 


10 


7.2E- 


10 


4.7E-10 


6.7E-08 


S140 








1.3E-07 






l.OE- 


10 


1.2E- 


10 


7.1E- 


11 






1.5E-10 


8.4E-08 




40 


20 


3.4E-07 






1.4E- 


10 


2.9E- 


10 


2.7E- 


10 






9.1E-10 


3.4E-08 



onto dust grains, recombination onto negatively charged 
grains is a negligible H;^ destruction process compared with 
this reaction. Therefore, for abundance in this regime 
we obtain: 



C/n 



(4) 



fciX(CO) 

where C, is the cosmic-ray ionization rate and n is the total 
gas density. Now combining Eq. (|3} and (|4]) we come to the 
following simple expression: 

C/n 



X(HCO"* 



(5) 



a(HCO+)Xe + fcgXg 

Estimates of the electron abundances based on this for- 
mula are presented in the last column of Table 1101 The 
cosmic-ray ionization rate was assumed to be e qual to 
C = 3 X 10"" s"\ a(HCO+) = 7.5 x 10~^ s'^cm^ jTurneil 
1 19951 ). The gas density was assumed to be n = 10 cm 
and [HCO+]/[H^='CO+] = 40. The fcgXg term leads to cor- 
rections less than 10% in the Xc values and we neglected 
it. Strictly speaking, in this way we derive more reliably 
the parameter nXc, i.e. the electron density, not abundance. 
However, for a better comparison with other results we will 
discuss further mainly the Xc values. 



5 DISCUSSION 

5.1 Variations of molecular abundances 

An inspection of the observational results and estimates of 
the physical parameters presented above leads to several im- 
portant conclusions: 

(1) There are strong variations of the relative inten- 
sities in the lines of different molecular tracers across the 
investigated high mass star forming regions. In Paper I we 
mentioned already striking differences between N2H+ and 
CS maps and the fact that the CS distribution follows that 
of the dust emission while N2H''' does not. Now we see that 
C^*0 and HCN, like CS, are good tracers of the dust emis- 
sion which presumably shows the total mass distribution. At 
the same time the behavior of HCO^ and HNC resembles 
that of N2H"'". Significant differences in distributions of var- 
ious species in par ticular sources have bee n noticed in some 
other studies (e.g. lUngerechts et al.lll997l ') but here we see 
systematic effects common for the sources in the sample. 



(2) There is no sign of CO and/or CS depletion in these 
objects (in contrast to cold dark clouds). The abundances 
derived for C^'^O from the J = 1 - and J = 2 - 1 tran- 
sitions more or less agree with each other. They are prac- 
tically constant within each source and among the whole 
sample (with small deviations which can be caused e.g. by 
temperature uncertainties) and are close to the "canonical" 
va lues derived earlier (e.g. X(C^^O) « 1.7 x 10~^ obtained 
bv lFrerking. Langer fc Wilsonlll98a ). 

(3) The N2H"'" and HNC abundances significantly de- 
crease with increasing ionization fraction (Fig. |8} . The best 
fit gives X(N2H+) (X X-i-^'^O-^ and X(HN^='C) oc 

with the correlation coefficients of \p\ ~ 0.85 for both de- 
pendences. At the same time the HCN abundance does not 
show significant variations. It i s interesting that the data for 
N2H"'' and Xc presented by iBergin et al. ()l999i ) are in ex- 



cellent agreement with the X(N2H"'") — X^ relation obtained 
here, as one can see in Fig.|8l where the open symbols are the 
X(N2H+) values derive d from the N 2 H+ an d C^*0 column 
densities presented by ^Be rgin et al.l (Il999ll. and assuming 
X(C^^O) = 1.7 X 10"^ ((Frerking et al.lll982l ). 

These correlations should not be affected much by the 
uncertainties in the abundances due to the LTE approxima- 
tion, given that both variables scale in a similar way if the 
non-LTE approach is adopted. 

(4) The derived ionization fraction in these objects in 
general increases towards the strong embedded IR sources 
which coincide with the main peaks of the dust and CS 
emission (by a factor of 2-3). The N2H+, HCO+ and HNC 
abundances decrease correspondingly in these areas. This 
is consistent with our results in Paper I where we found 
that the N2H^ abundance is systematically lower towards 
the CS/dust emission peaks which coincide with strong IR 
sources. However, it is important to emphasize that esti- 
mates of the electron abundance were performed under the 
assumption of a constant gas density. Our estimates of the 
gas density towards S255 IR and N give practically the same 
values. At the same time it is not unreasonable to expect 
density variations which can smooth the derived variations 
of the ionization fraction. 

It is worth noting that these variations of the ionization 
fraction refer to the values averaged over the line of sight and 
an increase of the electron abunda nce in vicinities of a young 
star can be significantly larger. lUngerechts et all (|l997 ) 
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Figure 8. Relative abundances of N2H+, HN^^C and H^^CN 
(from top to bottom) in dependence on the electron abundance. 
Our estimates are plotted by the fille d symbols. Th e open symbols 
correspond to the data presented bv lBergin et al.l |l99^. 



found similar variations of molecular abundances in the 
Orion molecular cloud and suggested that the abundances 
of molecular ions can be reduced by a higher electron abun- 
dance caused by UV radiation propagating in a clumpy pho- 
todissociation region. It is well known that in clumpy media 
UV photons can penetrate deep into de nse clouds leading t o 
significantly enhanced ionization (e.g. iBethell et aLll2007h . 
At the same tim e in their stu d y of t he ionization fraction 
in massive cores, iBergin et all (|l999l ) did not find any no- 
ticeable increase of the electron abundance from the edge to 
the centre of massive cores. This can be caused perhaps by 
insufficient luminosities of their sample sources. 

Thus, we believe that UV radiation from young massive 
stars in clumpy media can be the primary cause of the ob- 
served ionization enhancement. In ad dition. X-rays detected 
already from many massive stars (e.g. IZhekov fc PallalboOTi ) 
can be also responsible for this effect. 

(5) We see no dependence of the relative abundances 
on the velocity dispersion, as derived from the observed line 
widths. In principle such dependence could be expected in 
particular for N2 H^ if it really escapes perturbed reg ions as 
suggested e.g. bv lWomack. Ziurvs fc Wvckofj l|l992t ). 

(6) There may be a trend for decreasing N2H+ and in- 
creasing H^'^CN abundances with increasing temperature. 



However, this is based on only one point (W3 data) and 
cannot be considered reliable. 



5.2 Chemical implications 

Apparently the observed differences in molecular distribu- 
tions cannot be explained by molecular freeze-out as in low- 
mass cores. The kinetic (and dust) temperatures at the CS 
and N2II"'" peaks are similar in most cases and rather high, 
> 30 K, which probably makes freeze-out ineffective. This 
is confirmed by the absence of any indication of the CO and 
CS depletion as mentioned above. 

One possible explanation f or the observed chemical dif- 
ferentiations was proposed by iLintott et al.l |200^. They 
suggested that the enhancement of CS and reduction in 
N2II"'" abundance found in regions of high-mass star for- 
mation may be related to the high dynamical activity in 
these regions which could enhance the rate of collapse of 
cores above the free-fall rate. Consequently, high gas den- 
sities would be achieved before freeze-out had removed the 
molecules responsible for N2II"'" loss, while the high densities 
promote CS formation. 

However, this model has several drawbacks. In par- 
ticular, its predictions for some species (e.g. SO) are not 
supported by observations. Then, it predicts a decrease in 
the abundance accompanied by CS enhancement. Our data 
show no sign of the CS abundance variations. Nothing to 
say that an accelerated collapse itself has no more or less 
advanced physical basement. 

Our results presented in the previous section indicate 
that ionization fraction is probably more important in estab- 
lishing the steady-state N2H+ and HNC abundance in mas- 
sive cores. The only important process which forms N2H^ 



sive co res, ine only 
is (e.g. lTurneij|l995l ): 

H+-HN2 



N2H-^ + H2 

There are two main processes which destroy N2H^ 
reaction with the CO molecule 

N2H+ + C0^ HCO+ -I- N2 

and dissociative recombination 

N2H+ + e N2 + H 

N2H+ + e NH + N 



(6) 
the 

(7) 

(8) 
(9) 



which mainly lead to the form ation of N2 (90% of the total 
reaction), as recently found bv lMolek et al.l (|2007f ). 

Then, we have to take into account N2H+ recombina- 
tion onto negatively charged dust grains as in the case of 
HCO^. The steady state N2H^ abundance in this model 
will be given by the formula: 



X(N2H+) = 



fc2X(H+)X(N2) 



(10) 



k'iXiCO) + Q(N2H+)Xe + fcgXg 

where k2 is the rate of the reaction ([6|, fcs is the rate of 
the reaction ((7]) which is 8.8 x 10"^*^ cm'^s"'^ according to 
the UMIST database and q(N2H+) is the summary rate of 
the reactions ([8]), ((9]). Its value is significantly different in 
the UMIST and OSU databases. We accept the value of 
8 X 10~^ cm'^s"'^ at 30 K (E. Herbst, private communica- 
tion). This means that the dissociative recombination will 
dominate at Xc > 10"^X(CO), i.e. at > 10"'^ for the 
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standard CO abundance X(CO) ^ 10~*. This is close to an 
average electron abundances derived in this work. Therefore, 
in principle it is possible that the dissociative recombination 
of N2H+ can dominate at least in part of our objects. Then, 
if the Hg" and N2 abundances are more or less constant (as 
can be expected) we obtain for X(N2H+) the dependence 
on Xe similar to that presented in Fig. [8] 

It is less clear how to explain the behavior of the HCN 
and HNC abundances. Their formation pathways are very 
differe nt. HCN definitely forms from th e neutral-neutral pro- 
cesses (iTurner. Pirogov fc Minhlll997l ) 

N-h(CH2,CH3) ^ HCN (11) 

HNC forms via the distinctly independent sequence 

C+ + NH3 -> H2CN+ + e -> HNC (12) 

The main destruction processes for HNC are probably 
reactions with C^ and Hg" (|Turner et al.|[l997l ). Calculations 
of steady state abundances of these species require addi- 
tional chemical modeling. We note that the HNC formation 
is closely linked to NH3. Given that NH3 and N2H+ both 
form from N2, thus they are chemically related, one expects 
similar morphologies of HNC and N2H^, as in fact we ob- 
serve. 

5.3 Chemical indicators of massive protostars 

One of the most intriguing problems in studies of high mass 
star formation is the identification of massive protostars at 
the earliest phases of evolution. Some clumps from our sam- 
ple represent probably such protostellar objects. For exam- 
ple the SE clump in W3 area is associated with a water 
maser but has no embedded IR sources and/or UC H II 
regions. The HCO^ line profile shows a red-shifted self- 
absorption feature typical for a collapsing cloud. The mass 
of this clump from continuum data is about 70 Mq (Ta- 
ble O. Therefore, it can be a massive protostar on a rather 
early evolutionary stage. It is also very pronounced in N2H^ 
emission. Another example of this kind is the N2H^ emission 
peak in S187 (although we did not see signs of contraction on 
the line profiles). In the S255 area, the northern component 
which is dominant in N2H"'", is apparently much younger 
than S255 IR. These examples show that a relatively strong 
N2H"'" emission can be considered as an indicator of such ob- 
jects. The species which correlate with N2H^ (HCO"'' and 
HNC) can be also useful in this respect. 



6 CONCLUSIONS 

We presented and discussed here observations of five regions 
of active high mass star formation in various molecular lines 
(at 3 mm and at 1.3 mm) and in continuum at 1.3 mm. 
On the basis of these observations we estimated physical 
parameters of the sources and molecular abundances. The 
main results of this study can be summarized as follows: 

(i) The typical physical parameters for the sources in our 
sample are: kinetic temperature in the range ~ 30 — 50 K, 
masses from tens to hundreds solar masses, gas densities 
~ 10^ cm"'^, ionization fraction ~ 10^^. In most cases the 
ionization fraction slightly (a few times) increases towards 



the embedded YSOs. The observed clumps are close to grav- 
itational equilibrium. Our temperature estimates are sys- 
tematically lower (by a factor of about 1.5-2) compared to 
those obtained with NH3 observations. However, tempera- 
tures measured with CH3C2H are similar to dust tempera- 
tures, suggesting that the observed methyl acetylene tran- 
sitions better trace the dust than ammonia (1,1) and (2,2) 
lines. 

(ii) There are systematic differences in distributions of 
various molecules in regions of high mass star formation. The 
abundances of CO, CS and HCN are more or less constant 
and optically thin lines of rare isotopes of these species are 
good tracers of the dense gas distribution in these regions. 
There is no sign of CO and/or CS depletion as in cold low 
mass cores. 

(iii) At the same time, the abundances of the high density 
tracers HCO^, HNC and especially N2H"'", strongly vary in 
these objects. They anti-correlate with the ionization frac- 
tion (X(N2H+) oc X-i-3±o-3 and X(HN"C) oc X-« «±o-2) 
and as a result decrease towards the embedded YSOs. For 
N2H"'' this can be explained by dissociative recombination 
to be the dominant destroying process. This conclusion is 
more or less consistent with the data on chemical reaction 
rates. There is no correlation of these abundances with the 
line width. 

(iv) The described variations of the HCO"'', HNC and 
N2H"*' abundances make them potentially valuable indica- 
tors of massive protostars. In our sample there are some 
clumps which represent probably massive protostars at very 
early stages of evolution and they are very pronounced in 
the lines of these species, especially N2H"''. 
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